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SOME  METHODS  FOR  INCREASING  THE  RATI 

*  t 

OP  CALCULATION  OF  ELEXEBTAK3C  FUNC¬ 
TIONS  ON  DIGITAL  ELECTRONIC 
COMPUTERS 


Algorithms  for  table -polynomial  approximation  of  ele¬ 
mentary  functions  on  digital  ooaputers  are  dlsoussed.  At  the 
expense  of  a  certain  Increase  In  the  number  of  oomputer  memory 
cells  ocoupled  with  sub- routines  for  calculating  elementary 
functions,  these  methods  permit  a  significant  Increase  In  the 

j 

rate  of  calculation  of  these  functions. 


frj)  -  tt-  (.  2-  »*4  2 t  v 


-/ 


4, 


INTRODUCTION 


The  uaa  of  digital  electronic  computer*  aa  device*  to 
oontrol  complex  technical  lnatrumenta  leads  to  the  oreatloa  of 
computers  vith  a  high  operating  speed.  In  the  literature  auoh 
computers  are  called  oomputera  which  operate  In  the  real  time 
scale. 


Besides  the  Increase  in  the  rate  of  performing  individual 
operations  in  computers  which  operate  in  the  real  time  scale,  the 
problem  of  creating  mathematical  methods  which  permit  the  rapid 
performance  of  the  required  calculations  becomes  quite  real. 

Here  the  expediency  of  increasing  as  much  as  possible  the  rate 
ft?  such  baslo  elementary  functions  as  sinx,  coax,  tgx,  arcalnx, 
arctgx,  Inx,  ex  and^x  is  immediately  apparent.  For  many  devloes, 
where  the  digital  computer  Is  the  controlling  complex,  the  cal- 
eolation  of  elementary  functions  oan  occupy  up  to  40  or  even  60 
percent  of  the  time  required  for  the  solution  of  a  control  pro¬ 
blem.  Therefore,  with  a  reduction  in  the  time  for  computing  ele¬ 
mentary  functions,  the  total  solution  time  may  be  reduced  signi¬ 
ficantly,  and  the  freed  time  can  be  used  for  correcting  the  con¬ 
trol  of  the  aotual  Instrument. 

Algorithms  with  which  elementary  functions  are  calculated 
on  digital  computers  intended  for  the  solution  of  engineering  pro¬ 
blems  are  poorly  suited  for  computers  which  operate  in  the  real 
time  scale.  This  is  explained  by  the  fact  that,  for  usual  digital 


pTb-rr-  H-ltl  y/it 


-3- 


computers,  the  choice  of  algorithms  la  baaed  primarily  on  coa- 
alderatlona  of  the  apaoa  oocupled  bj  tbelr  programming  la  the 
memory  unit* 

Algorithms  whlob  are  presently  In  uae  are  oonstruoted 
In  sue*.  a  way  that  they  give  the  value  a  of  elementary  functions 
on  tbe  entire  Interval  of  variation  of  the  argument  with  a  given 
degree  of  aoouraey. 

Algorithms  used  for  tbe  oaloulatlon  of  elementary  funo- 
tlona  oan  be  divided  Into  four  groupst 

1)  Algorithms  which  permit  tbe  calculation  of  elementary 
funotlons  with  polynomials  which  approximate  the  given  function 

( polynomial  approximation) ; 

2)  Algorithms  which  permit  the  oaloulatlon  of  elementary 
functions  with  rational  funotlons  whloh  approximate  a  given  funo- 
tion  (rational  approximation); 

3)  Algorithms  whloh  permit  the  oaloulatlon  of  elementary 
functions  with  tables  whloh  contain  referenoe  values  of  the  funo- 
tlon  on  a  given  Interval  (table  approximation); 

4)  Algorithms  whlob  permit  the  oaloulatlon  of  elementary 
funotlons  with  the  method  of  calculation  of  values  of  the  funo- 
tlon  "digit  by  digit," 

The  enumerated  algorithms  for  calculating  elementary  func¬ 
tions  have  one  part  In  oommons  they  require  a  reduotion,  using 
appropriate  formulas,  of  the  entire  interval  of  variation  of  the 
argument  of  a  given  elementary  function  to  some  minimum  interval, 
where  tbe  approximation  Is  accomplished  by  a.  given  algorithm.  Ve 
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shall  consider  this  common  part  In  the  enumerated  algorithms  to 
he  identical,  and  we  shall  see  the  distinction  between  algorithms 
in  the  methods  through  which  the  calculation  of  the  funetlon  is 
accomplished  in  the  reduoed  Interval. 

In  widest  use  today  are  algorithms  whloh  permit  a  poly¬ 
nomial  approximation  of  elementary  functions.  Therefore,  these 
algorithms  will  be  of  greatest  interest  for  us. 

Algorithms  whloh  use  rational  approximation  are  not  yet 
in  wide  use,  since  in  the  majority  of  existing  computers  the 
operation  of  division  exoeeda  the  operation  of  multiplication 
in  the  time  required  for  its  performance  by  two  or  three  times. 
They  are  used  only  in  computers  in  which  the  operation  of  divi¬ 
sion  is  to  a  certain  extent  equivalent  to  the  operation  of  multi¬ 
plication  in  the  time  required  for  its  performance. 

Algorithms  which  use  tables  of  reference  values  for  cal¬ 
culating  elementary  functions  are  even  less  widely  used.  This  is 
-explained  by  the  faot  that  the  tables  are  constructed  to  use  the 
most  simple  interpolation  formulas  which  contain  a  minimum  num¬ 
ber  of  multiplications,  and  this  leads  to  extremely  large  tables 
of  reference  values  wben  the  condition  of  high  aoouraoy  of  cal¬ 
culations  is  imposed.  A  large  volume  of  tables  increases  the 
time  required  to  locate, the  necessary  values  and  requires  a  large 
memory  unit.  Tables  of  reference  values  of  functions  can  be  con¬ 
sidered  fully  adequate  for  lower  accuracy  of  calculation  not  ex¬ 
ceeding  0.5‘10'V  when  a  sufficiently  large  memory  unit  is 
available. 
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Algorithms  which  permit  calculation  of  alimentary  funo- 
tions  "digit  by  digit"  art  suitable  for  use  in  a  computer.  Hare 
certain  elementary  func tions,  for  example  lnx  and  •*,  can  be  eal- 
culated  in  the  time  of  one  multiplication.  The  inadequacy  of 
these  algorithms  is  tbs  fact  that  not  all  elementary  funotlons 
oan  be  calculated  with  this  method.  In  addition  each  elementary 
function  requires  expensive  additional  equipment  for  its  approxi¬ 
mation  by  this  method. 

Experience  In  the  operation  of  digital  eleotronlo  com- 

.  i  — 

puters  in  use  today  has  shown  that  algorithms  which  permit  poly¬ 
nomial  or  rational  approximation  of  elementary  functions  are  most 
suitable  in  practice. 

(Here  we  shall  examine  rational  approximation  in  the  re¬ 
duced  form,  l.o.  as  the  ratio  of  two  polynomials.) 

Thus  we  may  assume  that  in  computers  which  operate  in  the 
real  time  scale  these  algorithms  may  be  employed  with  certain 
variations  in  tbs  method  of  application.  With  an  Increase  in 
sub-routines,  they  may  be  redesigned  in  such  a  way  that  the 
approximation  of  any  elementary  function  in  the  reduced  Interval 
la  conducted  not  with  one  polynomial,  which  Insures  a  given 
accuracy  on  the  entire  interval,  but  by  a  certain. set  of  poly¬ 
nomials  applicable  on  sub-intervals  of  this  reduced  Interval. 

Here  the  degrees  of  the  approximating  polynomials  increase  from 
one  sub-interval  to  the  next,  and  only  in  the  last  sub-interval 
of  the  reduced  Interval  can  the  polynomial  have  the  fora  which  is 
now  used  over  the  entire  length  of  the  inverval  of  reduction. 
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Thus  a  combined  table -polynomial  algorithm,  as  It  were.  Is  being 
used  to  approximate  the  elementary  funotlons.  The  same  considera¬ 
tions  apply  to  rational  approximation  represented  as  the  ratio  of 
two  polynomials. 

I.  THE  TABLE-POLYNOMIAL  APPROXIMATION  METHOD 

The  possibility  and  expediency  of  employing  the  method 
of  table-polynomial  approximation  of  elementary  funotlons  may 
be  based  on  the  following  considerations. 

As  Is  generally  known,  a  computer  operates  with  digital 
Information  with  a  certain  degree  of  aocuraoy. 

This  acouraoy  Is  determined  not  only  by  the  error  In 
approximating  a  given  function,  but  also  by  the  faot  that  the 
computer  performs  the  neoessary  operations  on  numbers  which  lie 
In  a  fully  defined  range  and  are  discretely  distributed  within 
this  range.  The  accuracy  of  calculation  of  elementary  functions 
according  to  any  algorithm  Is  determined  also  by  this  discrete¬ 
ness.  For  values  of  functions  near  zero  In  modulus, polynomials 
which  give  a  good  approximation  near  the  values  of  these  func¬ 
tions  which  are  maximum  In  modulus  have  a  large  number  of  terms 
of  higher  orders  In  the  sub-intervals  near  the  zero  values.  How¬ 
ever,  the  terms  of  higher  orders  In  sub-intervals  near  the  zero 
values  of  the  function  take  on  values  which  are  not  caught  by  the 
computer^as  a  result  of  Its  limited  discrete  perception  of  these 
numbers,  but  .  require  a  large  number  of  multiplications. 
Herein  lies  the  basis  for  hope  that  the  entire  Interval  of  reduo- 
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tion  oan  Judiciously  be  divided  Into  smaller  sub-intervals  In 
whloh  an  approximation  of  functions  can  be  accomplished  with 
polynomials  of  lower  orders. 

Let  us  assume  that  we  have  succeeded  In  dividing  the 
Interval  of  reduction  (a,b)  Into  the  segments  A'p  Ag,  ...»  a^, 

suoh  that  At  H  A*  *  O,  \)  *  (a,  b). 

1  V  1*1 

Relative  to  the  time  required  for  their  accomplishment 
all  arlthmetlo  operations  performed  by  the  computer  can  be  ex¬ 
amined  from  the  point  of  view  of  the  time  required  for  the  per¬ 
formance  of  the  operations  of  addition  and  multiplication,  so 
that  the  operation  of  subtraction  Is  In  this  sense  equivalent 
to  the  operation  of  addition,  and  the  operation  of  division  oan 
be  thought  of  as  several  multiplications,  eto. 

Let  the  polynomial  which  gives  an  approximation  of  m 
particular  elementary  function  on  the  entire  Interval  (a,  b) 
require  m  additions  and  n  multiplications. 

The  time  expended  by  the  computer  to  accomplish  one 
or  another  algorithm  la  determined  as  the  number  of  arithmetic 
operations  of  addition  and  multiplication,  multiplied  by  the 
corresponding  time  required  for  the  performance  of  these  opera¬ 
tions  on  the  computer. 

Further,  let  nj  and  m^  be* the  number  of  multiplications 
and  additions  on  the  segment  A .  Here  we  assume  that  n^  <  «l4ir 
m^  4  *1+]/  Then  In  the  case  of  application  of  one  polynomial 
the  number  of  multiplications  on  (a,b)  is 
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/i 

lx{n)  -  n  ,  (I) 

* 

while  In  the  ease  of  application  of  several  polynomials  the  num¬ 
ber  of  multiplications  la 

x 

l2(n)  -  JZ  P(^i)nt,  .  (2) 

l "1 

where  P(*{)  are  the  probabilities  that  the  argument  for  whose 
values  the  elementary  function  is  being  calculated  will  fall  in  Aj, 

We  may  obtain  exaotly  the  same  expressions  for  the  mathe¬ 
matical  expectation  of  the  number  of  additions  when  calculating 
an  elementary  funotlon: 

i 

K.(m)  -  *  (3) 

K  i 

*aM  '  £  P(*i)-"l .  (4) 

The  time  required  for  the  accomplishment  of  the  algorithm 
on  the  oomputer  in  each  of  these  oases  will  bet 

i 

Tl  -Ti*i (n)  -  (fin  (5) 

T2  ^2*3  ^  '  (6) 

I 

I 

where  is  the  time  of  one  multiplication  on  the  computer, 

Tg  is  the  time  of  one  addition  on  the  computer. 

Let  ub  assume  that  a  polynomial  of  appropriate  degree 
which  approximates  the  given  elementary  function  on  a.  Is  given 
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by  lta  coefficients.  Then  the  selection  of  the  polynomial  reduces 
to  the  selection  of  the  appropriate  coefficients  from  the  table. 

Let  tbe  time  required  to  locate  these  ooefflolenta  on  the  tablo 
equal  .  Then  the  ratio  ^ 

1  .  ♦  «sV«J  *  **0  (7) 

<*in  ♦  *a«j 

I 

permits  an  evaluation  of  the  advantage  of  the  table-polynomial 
method  of  calculation  of  elementary  functions. 

An  evaluation  of  the  character  of  the  distribution  of 
probabilities  ?(x)  that  the  values  of  the  argument  x  will  lie 
In  the  corresponding  segment  of  the  reduced  Interval  Is  diffi¬ 
cult  even  In  the  case  of  control  units  of  relatively  slight  com¬ 
plexity.  We  shall  assume  that  B(x)  Is  determined  by  uniform  dis¬ 
tribution  and  that  the  oocurence  of  x  In  one  or  another  sub-interval 
of  the  reduoed  Interval  Is  an  independent  event. 

Then  from  (7)  we  easily  find  that 

I 

A  T°  .  (8) 

(b-a)  (nx^  *  nx%)  n-Cj  ♦  *t3 

I 

Thus  we  see  that  A  consists  of  the  two  parts 


(b-a)  (nx1  *  mxz) 
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The  ratio  Aj,  as  we  shall  show  subsequently  for  speoiflo 
elementary  functions,  can  be  made  leas  than  1  through  an  appro¬ 
priate  choice  of  approximating  polynomials.  It  will  be  near 
0.4  -  0.7. 

We  shall  oonslder  briefly  the  question  of  how  to  organise 
the  searoh  In  the  table  In  such  a  way  that  we  may  determine  the 
corresponding  A'  according  to  the  argument  x  reduced  to  the  baslo 
Interval  (more  precisely,  those  coefficients  a^  of  the  polynomial 
which  approximates  the  given  function  on  the  sub-interval  a- ) 

Z» 

with  the  condition  that  the  ratio  A^  -  be  as  small 

2  rWjj 

as  possible. 

We  may  consider  the  usual  program  methods  used  during  the 
choloe  of  one  of  k  possible  branches  of  the  calculation  process 
as  one  method  for  solving  this  problem. 

The  division  of  the  Interval  of  reduction  into  the  system 
of  sub-intervals  can  be  organized  In  such  a  way  that  the  sub¬ 
interval  A i  Is  determined  by  several  significant  digits  of  the 
binary  representation  of  the  argument  x,  or  that  the  argument  x 
is  associated  with  an  appropriately  seleoted  constant  prior  to 
the  selection  of  .the  sub-interval,  or  that  a  somewhat  more  com¬ 
plex  function  of  x  has  been  formed  which  contains  certain  charac¬ 
teristic  indications  for  each  sub-interval  A*  *  Then  the  choloe 
of  the  sub-interval  A,'  and,  consequently,  the  choice  of  an  appro¬ 
priate  approximating  polynomial  Pn  (x)  can  be  expediently  accom¬ 
plished  by  a  very  simple  deciphering  unit  which  will  control  the 
calculation  of  the  polynomial  with  the  required  number  of 
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coefficient*.  In  this  case  the  time  ta  required  for  the  search 
of  the  table  oan  be  made  very  email  in  comparison  to  the  time 
Tj  +  'tj;.  It  Is  true  that  In  this  case  certain  sub-routine*  will 
require  Individual  elementary  modifications*  which  do  not  com¬ 
plicate  the  design  of  the  compute r»ln  order  for  a  very  signifi¬ 
cant  advantage  In  rate  of  operation  to  be  obtained. 

In  essence*  we  may  provide  in  the  computer  a  special 
command  which  would  determine  the  appropriate  address  of  tha 
sub-routine  for  determining  Fn(x)  according  to  a  given  x.  Then 
this  command  could  seleot  several  digits  from  x  and  add  to  them 
a  certain  constant,  and  the  quantity  thus  derived  could  be  the 

address  of  the  sub-routine .  Ve  oould  provide  other*  more  complex 

*  % 

commands*  specially  adapted  to  this  problem*  whloh  would  permit 
a  maximum  reduction  of  the  time  required  to  search  the  Interval 
to  which  the  given  x  belongs.  This  would  permit  a  branching  to 
a  special  sub-routine  for  calculation  of  the  appropriate  approxi¬ 
mation  of  the  elementary  funotlon. 

II.  THE  TABLE-POLYNOMIAL  APPROXIMATION 
OP  ELEMENTARY  FUNCTIONS 

If  on  the  interval  £a,b3  we  are  given  an  arbitrary 
differentiable  function  f(x),  then  the  best  piecewise  approxi¬ 
mation  of  this  function  on  the  given  Interval  oan  be  performed 
with  polynomials  or  rational  fraction  functions  of  any  degree, 
with  the  assumption  that  the  Interval  £a,bj  la  divided  Into  a 
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series  of  smaller  Intervals. 

Vo  shall  perform  the  division  of  the  Interval. £a,bj  In 
one  of  the  following  two  ways. 

If  the  function  f(x)  has  the  property  of  evenness  or  odd¬ 
ness,  then  In  order  not  to  lose  these  properties  of  the  function 
we  shall  divide  the  Interval  |ja,b]  Into  a  series  of  nested  Inter¬ 
vals  with  the  common  point  a.  In  this  case  the  length  of  the 
nested  Intervals  Into  which  the  Interval  fa,bl  Is  divided  are 
determined  In  suoh  a  way  that  the  degree  of  the  polynomials  whloh 
approximate  the  funotion  f(x)  with  a  given  accuracy  Increases  with 
an  Increase  In  the  length  of  the  segment. 


If  the  funotion  f(x)  does  not  have  the  property  of  even¬ 
ness  or  oddness,  the  Interval  £a,bl  Is  divided  into  a  series  of 
small  non-overlapping  Intervals.  In  this  case  we  may  design  In 
advance  a  system  of  easily  calculated  best  polynomials  of  low  de¬ 
gree*  which  approximate  an  arbitrary  differentiable  function  on 
the  Interval  [x^,  xi4l?  *  ja,bj.  These  polynomials  will  be  re¬ 
quired  later  during  the  construction  of  formulas  for  approximating 
the  functions  /x  and  ez. 


We  shall  expand  the  given  function  f(x)  into  a  Taylor 
series  at  the  point  a^,  which  Is  the  mld-polnt  of  the  Interval 
fxl*  xlt£Ji 


f(x)-f(ai) 


(x-ai) 

+  ~zT 


f‘( <n, )  + 


(x-°l) 


(9) 
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Using  tbs  substitution  x-at  =  u  and  Cbebyohev 

polynomials  of  appropriate  degress,  and  returning  again  to  the 

'  » 

old  variable,  we  obtain  the  best  polynomial  of  the  necessary  de¬ 
gree  whl oh  approximates  the  funotion  f(x)  on  the  interval  £x^, 

Thus  the  best  polynoalal  which  approximates  the  funotion 
f(x)  on  the  Interval  jx^,  will  be* 


of  the  seoond  degree 


\ 


r(x)*ia  *  B.Jx-ai)  *  C2(x-ot)* 


(10) 


where  the  coefficients  of  the  polynomial  are  defined  by  the 
following  expressions* 


Az  “  Ni^a"  • 


Ea  -  [f(at)  --4sa>i*a  ), 

2*  2* 

c3  -  [—  x|V  *al. 

21  24 

*2 

t>2  "  C~3i f  "(a0  *~2?Sz)  **  ' 

sa  -JL  x*  ry(*0, 

41 


x3 


-  tn  -  *1 


-14- 


of  the  third  degree 


f(x)  <v  Ag  *  t'g'x -c^J  *  C8 (x-ai)2  *  ,  (11) 

where  the  eoefflolente  of  the  polynomial  ares 


A3  "  r>f(Qi)  *  f>S  18S3"  5,3  *3*3^ 


23 


b3  -  lf{at)  ♦ 1 ^±x3  j8), 

2*  24  3 

°»  ■’  l\>r<“0  -4<’-!s3  •  v. 

»» -  —  x|*3  1, 


-IV 


3 


3  "  I-,/  (at)  +-S9  Jx3  # 


28  3 


^3  4/^^  *4*3*3  ' 


23 


p  .  1  *  .vir,-  | 
"3  “giH  t  ' at '• 


S3  -  -l^l fn(*i)' 
x8  — 


m  ■  — ■  .r  r 
.  61 

*1+1  ~  xi 


of  the  fourth  degree 

f(x)  £  A4  *  B^x-at)*  Cifx-a^Difx-aiJt+e^fx-a^*  (12) 
where  the  coefficients  of  the  polynomial  ere: 
aa  “  lf(al}  -4?  *2  *4  * 


27 

B4  -  tf'(at)  *  Yex4S4  '■^4X4jr41' 

c4  -  tfj-  /"'(a*)  -i  >i  v  ft  x4*4 


i  [rlr^(ai)-L^84  h 

4  3!  23  23 
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*4  - 1^1*1)-  -p  A  h  -5-  f4  1  ' 

ir4  -  *  2*43*1  , 

#4  -  r^/7^; s4jx>  . 

^4.1X|  r^j, 

S4  -~  xtr713^  / 

61 

.  *1*1  *  xl 


\ 

I 

of  the  fifth  degree 

I 

f(x)  *.Ah+e5(x-ai)+CB(x-ai)2  *  Ds(x~ai)3*  es(x-at) *+  . 

*  T^x-arf*  ,  (13) 

where  the  coefficients  of  the  polynomial  are: 

tw 

A*  -  V<*0 -jA  s*  *Y* 

h  -if  (at)  -ix|  »8*JZ_jAr#j, 

28  2® 

c5-i^"W  ~  ^v|rxlv  . 
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1  it  *  a  3 
g6  't^T  f** (a0  --px6  S,  ♦  ; 

^5  't-r7  4XI  St  *2**»  ^  .* 

h  -  i—f^^O  *  2R5]  x|  , 

61 

H-^mo  , 

*tn  -  *i 

Xj  *  “'  • 

2 


We  now  proceed  with  the  direct  derivation  of  formulas 
for  specific  elementary  functions. 

Hereafter  it  is  assumed  that  all  formulas  constructed 
for  the  approximation  of  elementary  functions  Insure  the  derl~ 
vatlon  of  10  significant  decimal  places*  we  therefore  will  not 
evaluate  them  in  each  separate  oass. 

Here  we  shall  pause  to  analyse  the  formulas  given  in 
[,1^.  These  formulas  are  used  for  the  calculation  of  elementary 
functions  in  the  computers  BESM-1  and  BESM-II.  Some  of  these 
use  rational  approximation;  others  use  polynomial  approximation. 
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It  has  been  shown  In  that.  In  comparison  with  purs 
polynomial  approximation,  the  us#  of  these  formulas  permits  a 
significant  reduction  in  the  time  required  for  calculation  of 
elementary  functions.  Thus  a  further  refinement  of  these  for¬ 
mulas  is  desirable  for  the  purpose  of  obtaining  a  further  eoonomy 
in  calculating  time. 


The  following  considerations,  which  oonoern  the  calcu¬ 
lation  of  elementary  functions  in  the  smaller  sub-intervale  of 
the  reduced  Interval  with  polynomials  of  lower  degrees,  relate 
equally  to  any  other  formulas  which  use  polynomial  or  rational 
approximation. 

1.  Calculation  of  the  functions  slnx,  coax  and  tgx 


is  in  £lj,  the  expansion  of  the  function  ctgz  into  a 
Laurent  series  is  used  for  calculation  of  the  values  of  the  func¬ 
tions  slnx,  cosx  and  tgx: 

t 

zctgr  -  1  -  £  1  f3*-1  .  (14) 

K-l  (2x)  i 

which  converges  for  )z|  <  7f. 


Setting  a-  *  j,  we  obtain  the  expression 


I 

x 

K-l  (2 x) ! 

I 


(15) 


It  is  assumed  that,  before  calculation  of  the  tangent  of 
half  of  the  argument  according  to  formula  (15)»  the  original 
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argument  was  reduced  to  a  quantity  less  than  ^ ,  and  from  this 
quantity  a  whole  and  a  fractional  part  were  extracted.  The  ar¬ 
gument  x  in  this  formula  Is  the  fraotlonal  part  of  a  quantity 
less  than and  thus  varies  In  the  Interval  0(x$l. 

Ve  shall  divide  the  Interval  of  variation  of  the  lnde- 
pendent  variable  x  Into  four  partial  sub-intervals  with  the 
points  ^  ^  and  ^  . 

For  eaoh  of  the  sub-intervals  thus  obtained  we  oonatruot 
Its  best  polynomial  which  approximates  In  It  the  denominator  of 
the  function  defined  by  relatlon(15) .  Then  we  shall  have  the 
following  series  of  approximating  expressions. 

In  the  interval  £o;  ^Jthe  function  2tg-fka  approximated 
by  the  expressions 


where  the  coefficients  A^,  B^,  and  are  defined  as  follows! 

i 

-  [1  +  ~7M  Pl  ~  * 

*i 

-'TJ-S-  S&  *  4*1' 

h  -  ^ ;  A  . 

6.'  42 
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*1  *  IT  ~35  x? 

xi  • 

2* 


Formula  (16)  la  used  to  oaloulate  the  values  of  2t8$ 
for  anjr  x  on  the  Interval  Co»g  J« 

For  the  Interval  £o;  ^-Jwe  find  the  expression  which 
approximates  the  function  2tgX; 


8tg- 


2  A3  -  2  -  -  V®  ' 


(17) 


where  the  coefficients  Ag.  Bg,  Cg  and  Dg  are  defined  as  followst 


[1 


1  ,8 


'2  "  lA  2®  *3  *  J  , 

1  1  25  e  1  4 

5a  “  [iT  6*'  28  x3ff2  *  — 


1  1  25  5 

C  -  [ - +  —  x*  ^  — -Xs  *  ]  , 

2  41  30  26  2  8  23  3  » 


1  1  35 

D*  “  *-6T  42  “  24  ff3  *  ^3^ « 
115 

*2*^8!  16  +  T*2  Aa' 


_1_  _5  a 
*a  "  101  66  x3 

x3  *"*§“• 
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Formula  (,17)  la  used  to  calculate  the  values  of  2tg& 
only  for  the  argument  In  tha  interval  l£xt|,  alnoe  calculation 
of  the  values  of  the  funotlon  for  the  argument  x  on  the  Interval 
0*x$J-  la  aeooapll8hed  more  eeonomloall;  aQoordlng  to  formula  (16). 

In  the  interval  £0,|-  }tbe  following  expresalon,  whloh 
approximates  the  function  2tgj>  in  the  indicated  Interval,  holds s 


_  x 

2tgo 


2  ^8  “ 


~  D&fl  -  Xtffi 


(18) 


where  the  ooeffiolenta  A  ,  B,t  C  ,  D  and  £  are  defined  aa 

3  3  3  3  3 

follows l 


Aam  [l  xs°  V^5xi  'a1 


2® 


a,-  vs], 

3  2!  8  28  3  3  28  8  5 


1  1  105 


25 


"3 


■  [ - —  x®  Jf-  *  — -  xlifo)  / 

i  1  r\r\  /aft  ®  •  £6  “  * 


4i  30  28 


117.  35 

*  [ - - +  X3  “3  ”  Xg  Kg] 

3  6!  42  2s  3  24  8  3  * 

l  1  27  .  5 

*3  -  [— - 3  A  *3  *-z*»  J  » 

81  aO  23  3  8  2  s 

1  5  2 

JTo  -  ( — ; - 3 Ra)  xl  , 

3  101  68  8/  8 


- 


691 


3  121  2730  8  ' 

3 

‘!‘T- 


Formula  (18)  la  used  to  calculate  the  values  of  2tg& 
only  for  those  x  which  lie  in  the  interval! $x$jk 
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For  x  fro*  the  Interval  0$x<gihe  value*  of  the  funotlon 

,  »  ■  ■  .  • 

are  calculated  according  to  foraulaa  (16)  and  (17)*  :  - 


In  the  interval |b;  i]  the ln*±t^  expression,  which  approxi¬ 


mates  th*  funotlon  2tg  holds; 

I 

st *--■ 


(19) 


2  aa  -  bax*  -  cAx*  -  d4x6  -  b^x3  -  r4x^o 

I 

where  the  coefficients  cf  the  denominator  are  defined  as  follow*! 


1  l 

i.  -  ti  —  a4  +—xA)  , 

4  213  211  4 

1  1  49  9 

®4  '  C2lT  213*4  *  2®  *4  J* 

105 

R. - 1. 

4  2s  4 

*4  +  ] 

27 


t 

-  -  + 

41 

30 

1 

1 

l«T 

- - 

42 

1 

[  — 
81 

1  t  1 

30 

1 

r— 

10! 

5 

« 5 

1 

691 

12! 

2730 

1 

7 

7 

+  — Jf. 


11  a 

4  141  f. 


Formula  (19)  is  used,  not  for  the  entire  Interval  for 

which  It  was  derived,  but  only  for  values  of  the  argument  fro* 

the  Interval^.  $x$l. 

5 

Graphically  the  region  of  existence  of  each  of  expressions 
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(16),  (17),  (18)  and  (19)  and  the  Intervale  of  variation  of  the 
Independent  variable  for  eaeh  of  these  -expressions  oan  be  de¬ 
pleted  as  follows  t 


Value,  of  vgamm* 

Valua*  of  arguaMot  \ 
far  taywlw  (1S| 


Interval  of  de: 
of  cxptcMioa  (14) 


Interval  ofdefinkloa 
of  cxprcMloe  (17)  ; 


Value*  of 
far 


Value*  of  arguaeet 
far  epcptewlou  (U) 


Interval  of  defieitkja 
\  of  expreaalou  (U) 

,  Interval  of  definition  of  raprrjefan(lS) 


After  the  value  of  2tgjj.has  been  found,  the  absolute 
value  of  the  functions  slnx,  coax  and  tgx  are  easily  calculated 
according  to  the  foraulast 


,  *4 

slnx-  ■  . 

1  ♦!*»! 


cos x  * - 

1*  3L 
2 


(20) 


tgx 


sin  x 
cos  x 


/ 


For  calculation  of  the  functions  slnx,  cosx  and  tgx.  It 
would  have  been  possible  to  construct  still  other  approximating 
formulas  based  on  their  expansion  Into  a  Taylor  series  and 
Cbebyshev  polynomials  of  appropriate  degree. 

The  expansion  of  otgz  Into  a  Laurent  series  was  used 

.  ...  *  *  ,  t.j 

here  as  the  most  economical  for  the  simultaneous  determination 
of  all  three  functions.  Hence  we  have  demonstrated  a  method  for 
constructing  approximating  functions  which  accelerate  the  calou- 
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latlon  prooess  based  on  this  expansion* 

*  •• 

2 .  Calculation  of  the  function  Inn 

During  oaloulatlon  of  tbs  nstural  logarithm,  rtbor. 
argument  x.  Is  assused  to  be  given  in  normalised  fora*  l.s. 
In  tbs  fora 

x  -  2P*  xx  ,  (21) 

h 

where  p  Is  the  order  of  the  number  x,x£'ls  tip  Anginal  part  of 
the  nunber  x  whloh  Is  oontalned  In  the  Interval  jp^x  $!. 


Taking  the  logarithm  of  relation  (21)*  we  obtain 
lnx  -  pln2  ♦  lnxt,  (22) 

I 

where  lnx^  In  formula  (22)  oan  be  expressed  by  means  of  tbs 
series 


lnx. 


TV  * 1-1  *»♦* 


'r?4  xx  * 


2 


2n  +  l 


(23) 


I 

Using  relation  (23),  the  Identity 


lcXj 


J' 

Ini nxt  -  lni*t 


with  oertaln  assumptions  relative  to  the  particular  values  of 
fi ^ ,  and  Chebyshev  polynomials  of  appropriate  degree*  we  easily 
obtain  the  best  polynomials  whloh  approximate  the  function  lnx 
on  the  given  Interval  with  a  given  degree  of  aoouraoy. 


We  shall  divide  the  Interval  {] 


of  variation  of  the 
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Independent  Tarlable  x^  by  the  points  x^^?,  xj^and  X£  gtDtO 
four  partial  nested  sub- Intervale  with  the  common  point  x * ^  , 
and  In  eaob  of  the  intervale  thus  obtained,  including  the  entire 
lnt4fval  Itself,  we  uhall  construct  its  best  polynomial  whleh 
approximates  the  function  lnx. 

Then  for  the  Interval  Tlj  lllwa  will  have  the  following 

....  2  32 

relation  containing  the  beat  polynoalal  of  the  third  degree 

. ,  -v 

which  approximates  funotlon  lnx* 


lnx  -  (p-3)ln2  ♦‘glnl?  *  , 


(24) 


where  the  coefficients  are  defined  as  follows* 


AX  -  (2  fi  •>' 

B i  "  [3 *  2*  ' 


v-£*  4*^  * 

5  2 


VI 7 


u  -  ■ 


256  ♦  33YlT 
8 

arr  - 

M 

8 


Formula  (24)  Is  used  for  calculation  of  values  of  the 

function  lnx  for  all  x  lying  in  the  Interval 1.  $xy£ 

1  2  1  32 

The  following  relation,  which  contains  the  best  poly¬ 
nomial  of  the  fifth  degree,  is  obtained  for  the  lntervaljXj 
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In*  -  (p- 2)ln2  *  — ln5  ♦  iau  ♦  Ba  u*+  C#u*  >  (25) 

3 

vhara  tha  ooaffiolanta  art  daflnad  at  follows* 


4.  *  f  2  ~^-x!  ♦  —  x*  JTa  ) 
27  *  a* 


3  3 


'a*  **  '  2* 


2  3  4  7 

c2  ”  f - XT,  ♦  — jra  ) 

3  E  28  3  2*  3  ' 

2  1  3  9 

'a  -  f—  .-Xa  )  Xa  , 

V? 

Xa  "  lo  ♦  9f5  ‘ 

4 

"*  V  " 


Formula  (25)  la  used  for  calculation  of  tbe  values  of  the  funo< 
tlonlnx  only  for  those  values  of  x^  lying  In  the  Interval 

W  *■**  %■ 

For  z1$S’  It  la  more  expedient  to  use  formula  (24). 

1  32 

The  following  relation  with  the  test  polynomial  of  the 
seventh  degree  Is  constructed  for  the  Interval [1; 


where  the  ooeffiolents  of  the  polynomial  are  defined  aa  follows} 


12  *7  ? *  7** 

ft  s 

*•*  3  ¥ 


15  . 


2  7  27 

c3  -  f-  *  —xS  -  -  xf*3l  . 

5  2s  gi  - 
2  1  4  9 

*•*  ‘tt1*'*'8  ' 

2  1 

'a  '  (7*  Txl>  4  * 


X3  * 


9  2 

5|26 

104  ♦  2l'fi&  • 


8 


“s“VaT  ' 

n: 

1  o 


U  •- 


♦  SH 


Formula  (26)  Is  used  to  calculate  the  values  of  the 
funotlon  lnx  for  values  of  x^  lying  in  the  Interval^ 

The  following  expression  containing  the  best  polynomial 
of  the  ninth  degree  Is  constructed  on  the  Interval]]^;  1 

lux  -  (p-~)  ln2  ♦  A4u  *  B4u3+  C4us *  D4ur*  J4u®  •,  (27 


where  the  coefficients  of  the  polynomial  are  defined  as  follows: 

A4  m  £2  "£1  x\3*  , 

2  7  55  - 

B4  •  [~*  *  ~r§  J  , 

o  2*  2° 
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2  7  6  77  4  , 

4  5  2®  4  2®  4  * 

«  3  .  11  .  ,  ' 
*4  -  ff  ♦  -pjxj-  -jj-xjl*)  # 

2  5  4  U  , 

*4  '  • 

a  l 

Va  - 1 

.  - - 

4  V2  *  1 

iT  . 

X~¥L 

X*  2 

U  m  '  "  • 


Formula  (27)  la  used  for  calculation  of  the  valusm  of 
the  function  lnx  for  lying  In  the  lntervalg^x^l. 

The  regions  where  polynomials  (24)  -  (27)  approximate 
lnx  with  the  given  acouraoy  may  be  presented  graphically  am 
followsi 


Values  dt  argument 
foe  expression  (25) 

— — - \ 

Values  of  argument 
for  expression  (24) 


Interval  of  definition 
for  expression  (24)  / 

Interval  of  definition 
for  expression  (25) 


f 


Values  of  argument 
for  expression  (26) 


Values  of  argument 
for  expression  (27) 


Interval  of  definition 
far  expression  (27) 

Interval  of  definition 
\  for  egression  (26) 


From  the  very  construction  of  the  approximating  best 
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polynomial*  Introduced  hare,  it  follows  that  each  successive  poly¬ 
nomial  overlaps  the  preoeedlng  ones  and  permits  the  derivation  of 
the  value  lnx  with  the  given  acouraoy  in  those  intervals  mors  to 
the  left  than  the  region  of  variation  of  the  argument  for  whloh 
this  polynomial  was  oonstruoted.  Although  the  utilization  of  all 
of  these  polynomials  will  require  additional  loading  of  the  com¬ 
puter's  memory  unit,  it  will  nonetheless  lead  to  a  significant 
economy  of  maohlne  time  during  repeated  calculation  of  the  values 
of  lnx.  ^ 


3.  Calculation  of  the  funotlon  ez 

The  calculation  of  the  funotlon  e*  begins  with  the  re¬ 
presentation  of  this  funotlon  in  the  form: 


^  -  A  A1.  A’  - 


(28) 


where  f  1  is  the  whole  part  of  the  index  X  ,  and  {_iJls 
ln2  ln2  152* 


the  fractional  part  of  the  index 


ln2 


The  part  of  the  quantity  ex  which  is  expressed  by  the 
r_jt  -l 

faotor  we  assume  to  be  calculated  and  ready  to  add  to  the 


other  part  of  e  ,  l.e.  to  2 


,  which  remains  to  be  calculated. 


Where 


Here  we  note  the  applicability  of  the  equality 
2  -  •  m  •  '  (29) 

«-{iTT)ln2- 

Clearly  the  variable  V  varies  within  the  limits  JvJ^ln2, 
so  that  the  quantity  {  ^varies  between  0  and  1. 
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The  Halts  of  variation  of  the  index  funotion  oan  be  re¬ 
duced  further  If  we  multiply  and  divide  it  by  VF.  Then  we  obtains 


In* 


12 


(30) 


where  u  .  {-^l^  -  ^  varies  within  the  limits  |u|  <  1|£.  * 

Calculating  eu,  we  at  the  tame  time  find  the  value  of  tbe 
quantity  2-(ld^}  according  to  (29)* 

The  values  of  the  funotion  eu  can  be  calculated  in  two 

way*. 

First  method. 

As  the  Initial  aeries  for  calculation  of  the  funotion 
eu  we  take  its  expansion  Into  a  Maelaurln  series: 


U  li*  it  K*  tt® 

*a  '~5'" 

We  shall  divide  the  Interval  (p;  ^5^3 


(31) 


of  variation  of  the 


modulus  of  the  Independent  variable  u  by  the  three  points 

321b2  and  into  four  nested  Intervals  with  the  common  point 

1024  24 

0,  and  we  shall  construct  for  each  of  these  Intervals  tbe  best 

polynomial  which  approximates  the  funotion  eu. 

.  / 

Using  series  (31)  and  Chebyshev  polynomials,  for  the  in¬ 
terval  [0;  §|§3*e  obtain  the  following  best: polynomial: 


eu  -  ij  ♦  B^u  >  CjU8  , 


(32) 


where  the  coefficients  are  defined  as  follows: 

si.  -  ti  — * — ■*:  V 
2*  2s 
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°1  ■  4'  7  ‘Jv  "1 1  • 


h  -  '"a,  >i  . 

r  -  *  * 
s,  -  —  Ju 
1  41 


sl 

6l"  * 

3ln2 

**  “  4096  ‘ 


The  beat  polynomial  (32)  la  used  for  calculation  of  the 
▼alues  of  the  funotlon  eu  for  all  values  of  u  In  the  Interval 

H8R1 

The  following  beat  polynoalal  la  constructed  for  approxl- 
nation  of  the  funotlon  eu  on  the  Interval  [b; 

•*  ea  -  Aa  *  Bzu  *  Caua  *  Dqu*  *  ,  -  l®®' 

where  the  coefficients  are  defined  as  follows: 

1  6  l 

ka  “  ^  ~~$r  xzsz*  • 

Ba  “  t1  xPa  "*^4X3  *3  ^  >  , 

1  1  4  9  a 

Ca~  &  *  l2  NjV^a**1  * 

-  r  1  1  3-  5  - 

D 2  *  [*3i  xa*a  *-$*2  1  ' 

*s‘  'T'  "P*8  S>‘F,J  ’• 


31 


*  (—*  “A)  ** 

5?  2*  .  * 

•  l 


B  1  * 
“Ti^8  ' 

1  2 

*"“5*  * 


X3 


391n2 

1024 


The  polynomial  (33)  la  used  for  calculation  of  the  val- 

* 

ues  of  eu  only  for  those  values  of  the  argument  whloh  lie  In  the 

Interval  I 

3ln2  ^  39ln2 

-  <  “  <  -  ♦ 

4096  1024 

t  ’ 

For  the  Interval  [b;  1  obtain  the  following  best 

polynomial  to  approximate  the  function  eus 


8“  '  V  V  +  V*  *  D9U*  *  V*  *  V5  *G3U*:  (34) 


where  the  ooefflolenta  are  found  from  the  relations: 


i3  “  ^  *  2»>,f  S3~  9.7 


Z7 


C*-[T§  **  ***»>' 

1  8  1.  7  ft 

^3  t-  *  - 1*3*8^  » 

31  •  26  3  2s 
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1  26  .  8  , 

'•  '  <1  *2»  *»*•  'H  lj,»  1 

'*■  firslS''**  ^“*! ' 

1  35  . 

1  1  > 

'»  '  '~\  *•  ' 

*|v>l  , 

*•  -i*  ■ 

1  * 

So  *  X«  / 

3  101  3 

51a  2 


The  polynomial  (34}  la  used  for  oaloulatlon  of  tha  val¬ 
ues  of  eu  for  values  of  the  argument  eontalned  In  the  Interval 

1b2 

For  the  Interval  [b;  -jp']  we  obtain  the  following  best 
polynomial  to  give  the  values  of  the  function  eu j 


eu  «  X4  *  B^u  +  C  4  u3  ♦  D4 j*3  ♦  *  P+u  3*  ♦ 

♦  ♦  «/4u®  ,  (35)  ; 

where  the  ooftffloldnta.of  the  polynomial  are  defined  byi 
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i4  -  tl  -  *  g  *4*4  1  » 

1  1  8  9  . 

.1  9  •  25  - 

°4  ■  C3?  *  Ja*s*  "15^  f4  1  ' 

1  5  .  15  4 

**  "  ai"  ’*?  ' 


*4* 


1  105  25.  4 

[ - x?s4  ♦ —  v;  i4  ], 

4!  28  4  4  26  4  4 


1  7  4  27 

^Te^'F^41' 


B.  - 


11  9 

X4?4  *  "p  *4  ^  ' 


I.  - 


- 


B.  - 


1  27  -  5 

l— - -  x4  S4  ♦  —  tA 

8J  g3  44  2 

l  1  „  .  a 

li  *  ?V  l<  '  . 
^  ^  l*  ' 
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1  a 

12!"  V 
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Polynomial  (35)  la  used  for  oaloulatlon  of  the  values 
of  au  only  for  those  values  of  the  argument  lying  In  the  Interval 


51n2  ^  _  ln2 

-  $  u  < - - 

24  2 

• ! 

Seoond  aethod.  ! 

He  shall  divide  the  positive  part  of  the  Interval  of 
variation  of  the  Independent  variable  u  Into  8  equal  parte  and 
determine  the  midpoint  of  eaoh  of  these  segments. 


In2 

and  that  ^  s  32  identically  for  all  of  the  segments,  ve  shall 
find  the  values  of  the  funotlon  eu  and  Its  successive  derivatives 
at  the  points  a^  (36). 

Through  estimations  ve  find  that,  with  an  aoouraoy  not 
less  than  0.3*10  a  polynomial  of  tbs  fourth  degree  approxi¬ 
mates  the  function  eu  on  eaoh  of  the  segments. 

Thus  to  find  the  values  of  the  funotlon  eu  on  eaoh  of  the 
eight  segments  ve  oan  use  a  best  polynomial  of  the  form  (12))flnd- 
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lng  the  coefficients  of  these  polynomials  according  to  the 
same  formulas  as  for\(12),  Bara  we  naad  to  find  only  the  Taluaa 
of  tha  funotlon  and  lta  derivatives  to  the  eighth  order,  Inclu¬ 
sively,  at  eaoh  of  tha  polnta  a^  (1  ■  1,  2,  ...,8). 

With  thla  aathod  of  calculating  au  5  addltlona  and  4 
aultlplloatlona  must  ba  expanded  In  tha  derivation  of  tha  Taluaa 
of  the  funotlon.  Tha  neoeesary  Interval,  depending  on  tha  value 
of  u,  la  found  by  means  of  3  ooaparlaona. 

Tha  foregoing  discussion  relative  to  u>0  ralataa  equally 
to  u<Oy  which  also  require  approximating  polynoalala  of  tha  4th 
degree.  Thus  one  more  conditional' step  la  added  for  recognition 
of  the  algn  of  u. 

Thus  96  constants  of  16  approximating  polynoalala  must 

e 

be  stored  In  the  computer's  memory  unit  to  determine  the  values 
of  eu  for  any  u. 

After  the  quantity  eu  has  been  determined,  the  value  of 
the  funotlon  ex  la  found  by  means  of  relation  (26). 

4.  Calculation  of  the  square  root  Vx 

In  the  calculation  of  the  square  root  It  la  assumed 
that  the  argument  la  always  given  In  normalised  form,  l.e.  that 
aljl  values  of  x^  lie  In  the  Interval  |-^x^<l.  Then  calculation 
of  the  square  root  of  the  number  x  reduoes  to  the  extraotlon  of 

the  square  root,  of  the  two  faotorsi 

*  * 

*  . 

»  f 

i 

Extraction  of  the  square  root  of  the  first  faotor  poses 
no  particular  problems  slnoe 
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w 


2% 

2v*l. 


To  extract  the  square  root  of  the  second  faotor  we 
divide  the  region  of  definition  of  the  argument  into  16  equal 
segments, and  oonstruet  for  eaob  of  them  the  best  polynomial  to 
approximate  the  function  f(x)  =  Vx  in  that  interval* 

We  shall  denote  bj  a^  the  midpoints  of  the  segments* 
Then,  with  the  above  division  of  the  total  interval*  we  shall 
have  the  following  values  for  the  points  a^t  . 


33 

41 

49 

57 

ai  "  64  * 

“»  ’"54  ' 

°1S  "64* 

35 

43 

51 

59 

a  64 

**  'I?' 

9 

°14  *  54 

37 

45 

53 

61 

fl.  rm  — 

u3  " -  * 

64 

a»  - - / 

54 

aii  t 

64 

16 

64 

39 

a*  - - , 

4  64 

47 

08  'IT' 

55 

Ol  f*x  ‘  7 

64 

63 

o16»  1 

64 

The  distribution  of  the  Intervals  in  which  polynomials 
to  approximate  the  function  will  be  constructed  can  be  deploted 
graphically  as  follows: 

,j(. 

I  _9.  J9.  J5.  _21_  _n_  _23_  _3_  _25_  J3_  2?  7  29  15  31  4 
T  32  U  32  8  32  14  32  4  32  16  32  T  ~W  "IT  1 

*:  Bi  °3  *4  °»  “•  Br  *»  “»  “10  *u  #i»  “i*  *14  °i§ 


In  accomplishing  this  division  we  at  the  same  time  have 
determined  the  quantity  A  =  ~  -1 .  in  the  present  oase  A 

is  the  same  for  all  segments  and  equals  • 
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Calculating  the  valuta  of  the  derivatives  for  tha  points 
a^  (1  *  1,  2,  . ..,  16)  froa  (37)  according  to  the  femulaat 


r(at)  , 

„  1  * 
f*(a, )  «- —  - — 

l'  2  \pi 


1  1 


2s  (V^73'> 


fntct;  ■  -2i 


2® 

111!  1 


i  f™(ai)  m  - 


tfF2' 


t^H) 
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1 

1  2*  (Vop*  * 

SI! 

1 

1  “  2* 

711 

i 

2® 

131! 

■ - „ 

i 

2® 

15!! 

i 

2®  bMOvr  * 


(36) 


and  oatlmatlng  tha  size  of  the  remainder  tame,  wa  find  that, 
vlth  an  aoouraoj  not  lean  than  0.35*10“12,  the  funotlon  Vx  la 
approximated  In  the  flrat  fire  segments  froa  the  left  with  hast 
polynomials  of  the  fifth. degree  of  the  fora  (13)*  The  coeffi¬ 
cients  are  determined  In  the  same  way  as  for  the  polynomial  (13)* 
Here  the  values  of  the  function  and  its  derivatives  are  deter¬ 
mined  aooordlng  to  formula  (38)  for  eaoh  point  a^(l  a  1,2, 3*4, 5) 
from  (37). 

With  an  acouraoy  not  less  than  0.8*10~^,  tha  funotlon 

r 

_ .  t 

f(x)  *  yx  Is  approxlmatad  on  all  of  the  remaining  segments  by 
4-th  degree  polynomials  of  the  fora  (12),  where  the  ooefflolente 
are  determined  In  the  same  way  as  for  polynomial  (12).  Tha 
values  of  the  derivatives  are  found  aooordlng  to  formula  (38) 
for  eaoh  point  a^  (1  =  6,  7,  ...»  16)  froa  (37). 
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Thus,  In  dividing  th«  interval  C|*  lj  Into  16  partial 
segments,  we  ara  abla  to  oaloulata  tba  root  of  x^  with  ona  of 
16  polynomials,  requiring  In  tha  worst  case  6  additions  and  5 
multiplications,  and  In  tba  bast  ossa  5  additions  and  4  multi¬ 
plication*.  Finding  tbs  necessary  Interval  raqulrss  4  compari¬ 
sons. 

In  tba  oosputar'a  memory  unit  101  ooaff lolants  mu at  bs 

stored . 

It  la  posalbl#  that  the  following  method  of  calculating 
tba  values  of  tba  square  root  with  an  aocuracy  to  tan  dsolmal 
places  would  be  more  economical. 

Aa  before,  we  divide  tba  interval  of  definition  of  tba 
independent  variable  Into  16  consents. 

We  shall  us*  second-degree  polynomials  of  the  form  (10). 
Sstlmates  Indicate  that  tba  accuracy  of  the  square  root  calcula¬ 
tion  with  these  polynomials  will  not  be  lower  than  0.3*10"^  on 
any  of  the  16  segments.  Adding  one  Iteration,  which  doubles  tba 
accuracy,  wa  obtain  the  value  of  the  square  root  in  4  additions, 
3  multiplications  and  1  division. 

If  calculation  of  the  value  of  the  square  root  Is  per¬ 
formed  with  second-degree  polynomials  with  one  Iteration  for  the 
8  left  segments,  and  with  fourth-degree  polynomials  for  the  8 
right  segments,  the  number  of  coefficients  to  be  stored  reduoas 
to  60. 

If  the  value  of  tha  square  root  In  all  16  Intervals  la 
calculated  with  second-degree  polynomials  with  ona  subsequent 
Iteration,  the  number  of  coefficients  to  be  stored  will  be  64. 
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5.  Calculation  of  the  function  arotgx 


The  function  arotgx  can  be  calculated  for  all  number* 
x  which  can  be  represented  in  the  given  ooaputer. 

However,  the  foraulas  with  whloh  arotgx  la  oaloulated 
bold  only  for  |x| 4  1.  Therefor,  if  |x|  >1,  calculation  of  tha 
function  arotgx  is  accomplished  with  the  formula 


i 

arctg*  -  -  arct£— .  (39) 

1 

l.e.  it  is  reduced  to  the  calculation  of  the  arctangent  of  an 
argument  leas  than  unity. 

The  point  A  Isolates  from  the  Interval  -l$x<l  the 
Interval  |x|£Xsuch  that  the  polynomial  obtained  by  expanding 
the  funotlon  arotgx  into  a  Maolaurin  series  and  rotated  by  means 
of  Chebyahev  polynomials  gives  the  values  of  the  arctangent  on 
the  Interval  |xl  i  X  with  the  required  accuracy  when  the  polyno¬ 
mial  of  a  given  degree  is  used.  Thus  we  are  aile  to  calculate 
the  values  of  the  arctangent  on  the  Interval  lx)*  X.  To  calcu¬ 
late  arotgx  on  the  Interval  fA;  1]  ,  we  use  the  familiar  formula 
for  the  sum  of  two  arctangents: 


U+Vf 

arct gu  *  arctj?ot  -  aret^y— -  • 


(40) 


If  we  use  the  substitution 


1*1  - 

H*l«t  ' 

from  relation  (40)  we  easily  obtain: 

.  I. 

arctgu  +  arctgVj  -  arctjfx . 


(41) 


(42) 


Thus,  to  oaloulate  th«  values  of  arotgx  at  any  point  of  tha 
lntonral  *ilx)£l,  for  a  given  x  In  this  Interval  ve  suet  eal- 
oulate,  according  to  foraula  (41),  the  u  vhloh  corresponds  to 
It,  taking  the  value  of  v^  neoeesary  for  tbla,  and  substitute 
tbla  value  of  u  In  plaee  of  x  in  the  derived  beet  polynomial. 
Then  ve  add  the  value  of  arotgu  tbua  obtained  with  tbe  value 
of  aretgv^  calculated  earlier,  vhloh  agreee  with  relation  (42). 

In  order  to  be  able  to  use  tbe  beat  polynomial  vhloh 
holds  on  tbe  Interval  |x/<Ato  oaloulate  tbe  values  of  tbe 
function  on  the  Interval  A  <  fx|  £  1,  we  must  Insure  that  tbs  u 
vhloh  corresponds  to  a  given  x  Is  not  greater  than  A.  Hence 
tbe  Interval f  A;  l"]  la  divided  Into  the  k  Intervals  £x^;  X2I, 
[x2*x3l*  •••»  0V*7.  ln  oaob  of  which  |u)U<  The  dividing 

points  are  date mined  from  tbe  foraula 

1 

x{  *  Vi 

"  1  .  j.  V  /  #  (43) 

Tl 

where  x^  ■  A»  In  eaoh  of  the  Intervals  (x^;  x^^J  the 

constants  v^  are  found  from  the? foraula 


Vj  »  — . —  • 

1  \-xiXx 

i 


(44) 


The  values  of  arotgv^  are  oaloulated  In  advanee. 

Now  we  shall  consider  the  construction  of  approximating 

polynomials  on  the  Interval  )x|  <  A. 

V3T1 

If  we  set  A^  *  jTi*  the  best  polynomial  vhloh  provides 
the  Indicated  aoouraoy  on  the  Interval  Ixl^A^  will  be 


aretgx  -  ixx  -  Sxx3  *  Cxxs  -  Dxx7  +  Exx9  -  FjX11  , (45) 
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where 


n  1  .14  18  1*  „  . 

1  “iu*1  .  ^ 


-  [--—-4*  ♦  —  xi%], 

1  3  210  3  n  210  1  ■  * 

1  1  63  10  91  „ 

ci  -  i - -  xi°-  — xflxl  , 

5  2®  5  2 

n  4  15  a  ^  39  •  , 

7  27  2* 

_  J,  1  55  .  65  4 

1  '  fsT"?'3  ^  “24  X1  *ll  I 

1  3  13  _ 

TThh* 


2* 


The  Interval  >  £  )xj  &■  1  le  not  divided  Into  smaller  Intervale, 

For  It  ri  -  ^3  a*»d  arotsvj  *  ^  • 

The  degree  of  the  polynomial  whloh  approximates  tbs 
function  nrctgx  on  the  Interval  lx) £  x.oan  be  made  lower  than 
the  degree  of  polynomial  ( 45)  If  we  set  ^  .  Then  the  best 

polynomial  which  approximates  arotgx  on  the  Interval  lx!  $ 
will  be 

arctgx  -  AgX  -  B^x*  *  CgX*  -  D^e7  ,  (45) 


where 


1  1 A  9  O 

J*~  U~#b  **  ‘ze  VV  '  . 

15,  15  , 

*»  ”  [r_  ■  i*  v*J » 

1  7  ,  27 

c*“  "  -gi  • 

_  ,*  1  *  9  2 

a  *  C7 -  &  Xa  “"P  XslaJ  * 

1  1  s 

**  *  7’IF**  * 

For  calculation  of  arotgx  on  the  remaining  Interval, 

the  Interval  £  A-gjl }  la  divided  Into  the  three  Intervals  (pigjx^J , 

fxgixjl*  r*3 » •  The  values  of  v^,  Vg  and  v^  are  found  from 

formula  (44).  Then  the  values  of  arctgv1#  arotgVg  and  arotgv^ 

are  oaloulated.  The  values  of  arctgx  on  the  Interval  fAgjll  are 

calculated  with  formulas  (46),  (40)  and  (41)  using  one  of  the 

three  pairs  of  final  constants. 

Without  sacrificing  accuracy,  ve  may  lower  the  degree  of 

the  approximating  polynomial  still  further  If  we  take  A,  •  ^7  . 

v  lo 

In  this  oase  the  best  polynomial  which  approximates  the  func¬ 
tion  arctgx  on  the  Interval  fxf^A^  with  the  same  aocuraoy  will 
have  the  form 

where 


I,,.- 

arctgx  -  Agt  -  flgxS  ♦  Cjc&  , 


(47) 


1  .8  7  « 


A3  “  ~  jgS  X3  “  96  *3  *3  ^ 


2« 


1  1  5  .  7  4 

E*  "  l~~  7s  "7  %3  ' 

C  26  3  23 
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-*-4 
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1  1 

’77 


7  I, 


To  oaloulate  the  function  on  the  remaining  Interval,  we  divide 
the  segment  fl^all  Into  6  segue nt si  £ Xy,x^J,  (x^ixjjt  ...» 

(x^;  ll»  for  whloh  the  values  v1§  v2,  . ...  v^  and  arotgv^, 
arotgv2,  ...»  arctgvg  are  found.  When  these  values  are  used, 
the  value  of  u  does  not  exceed  X^  on  eaoh  of  the  segments,  and 
formula  (47),  together  with  formula  (40),  will  give  the  values 
of  arotgx  on  the  entire  Interval 

A  polynomial  of  lower  degree  will  approximate  the  func¬ 
tion  arotgx  on  the  Interval  |x|<Avlth  the  same  aoouraoy  if  we 
se t  this  oase  the  best  polynomial  whloh  gives  the 

values  of  the  funotlon  arotgx  on  the  Interval  JxfO^  will  bet 


where 


arctfjx  -  , 


(48) 


4  “  1  “  2«X*  "  3' 

;  1  .  5  2 

h“  £3  X4  '  23  x4*43' 

1  1  8 

*4  -7-?  ^4  * 

Thd  Interval  [  X^; 1}  is  divided  bj  the  points  x^,  determined 
by  formula  (43),  Into  the  29  segments  [X^ix^J,  fXgjXj}, 

for  whloh  the  values  of  Vj,  v2,  ...»  Vg^  and  arotgv^, 
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arctgv2,  ••••  *rotgv2 are  oaloulated.  With  these  oonstante  we 
■ay  mabe  the  absolute  value  of  u  leas  than  or  equal  to  X^  oa  ahy 
of  the  aegmenta,  and  compute  the  values  of  arotgx  on  the  entire 
Interval  [0;1'J  using  fornulas  <49)  and*  (40). 

Finally,  the  degree  of  the  approzlaatlng  polynomial  oan 
be  reduoed  to  unity  If  we  reduoe  further  the  Interval  on  whloh 
this  polynonlal  la  defined. 

Setting  Xg  «  gjjg  »  for  any  x  on  the  Interval  |x) 4  Xg 
we  obtain  the  values  of  the  funotlon  arotgx,  accurate  to  ten 
deolaal  places,  using  the  following  best  polynonlal  of  the  first 
degrees 

I 

Atctgx  -  Ajc  ,  (49) 

/  1 

1 

I 

where 


1  .  5  a  3  2.1  1  . 

*6  "  ^  ~  ~  -^4X6^6^'‘’ "«2X5^j  "*  ,,3X6 

-  ,1  1  v3  .  8 

* 6  m  ( - ~5X6  *  X6  • 

5  2s 


g2  *8^  • 


Clearly,  It  Is  not  posslDle  to  obtain  a  more  simple 
polynomial  than  polynomial  (49)  for  calculation  of  arotgx  at 
any  point  of  the  interval  |x|£X« 

The  Interval  £  1^; l^  1b  divided  by  the  points  x^  Into 
more  than  100  segments. 

If  for  these  segments  we  determine  the  constants  v^ 
and  arotgv^  from  formula  (44),  then  we  may  use  formulas  (49) 
and' (44)  to  calculate  the'  values  of  arotgx  on  the  entire  Inter¬ 
val  [6;ll. 
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6.  Calculation  of  the  function  aroalnx 


During  calculation  of  tbc  funotlon  aroalnx.  It  la  aaauncd 
that  tba  argument  la  contained  In  tbc  Interval  -If  xt  1,  l.e. 
that  only  tbc  principal  values  of  tba  funotlon  arc  calculated. 

Aa  the  Initial  8arlea  we  take  the  aerlea 

•  '  (2n— 1)1 1  jc3"*1  . 

arcsicx  -  *  j^li  ~^T  ’  (*») 

I 

defined  on  the  Interval  )x|  $1. 

We  shall  divide  the  Interval  of  variation  of  the  inde¬ 
pendent  variable  Into  four  partial  sub-intervals  by  the  polnta 
1/10,  2/5,  and  3/^,  and  construct  In  each  of  these  sub-inter¬ 
vals  [b;l/10l,  (jb ; 2/5l ,  [P;3/4J  and  £o;l"]  the  best  poljnoelal 
which  approximates  the  funotlon  aroalnx  In  It. 

Using  series  (50)  and  Chebysbev  polynomials  of  appro¬ 
priate  degree,  for  the  Interval  £0;l/l03  we  obtain  the  follow¬ 
ing  best  polynomial: 

arcsinx  -  ♦  3jX3  *  *  0^x7  ,  (51) 


where  the  coefficients  are  defined  as  follows: 

Ai  “  +-^0x6i*i  • 

r1  5  4  15  4  , 

Bi  -  ^  Hh  * 
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The  beet  polynoalal  (51)  le  used  to  o&leulate  the  valuea 
of  the  funotlon  aroslnz  on  the  Interval  [bjl/lO}.  For  the  inter- 
val  £o ; 2/51  ,  a  test  polynomial  of  higher  degree  la  oonatruotedi 

a  I 

arcs  in*  »  A^c  *  *  D^c*  * 

V  *  V11  (52) 

where  the  ooefflolenta  are  defined  aa  follows: 


ljra)  • 


1  12  15 

A2  "  ^  gieX3  S2  *  gl4 

1  3  35  _ 

Bz  *  f^'+'paxl0,a  --£n>xl°M2  1  »  . 

_  f3ll  1  21  8  189  8  , 
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5H  1  33  235 
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Polynomial  (52)  la  used  to  oaloulata  raluaa  of  the 
function  aroalnx  not  or  the  entire  interval  of  definition,  but 
only  for  those  values  of  the  argument  x  which  lie  in  the  Inter¬ 
val  1/10  £■  lx  I  £  2/5.  For  the  values  |x)  <  1/10,  calculation  of 
the  function  aroalnx  with  the  aame  aocuraey  la  muob  more  eaally 
accomplished  with  polynomial  (51)* 

A  best  polynomial  of  21at  degree  la  oonatruoted  to  ap- 
proximate  the  function  aroalnx  on  the  Interval  C0;3/4l  * 

arcsinx  -  K^x  *  Ijc3  ♦  t^x6  *  K^x7  * 

♦  r^9  ♦  fgx11  ♦  {jx1*  *  fgxic  ♦  f  l .  (63) 

*  F10*M  *  lu*» 

where  the  coefflolenta  are  defined  by  the  expreaalonst 

‘  1  23 

2d  20 
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J  <|x|  <  1,  we  oould  oonatruot  the  beet  polynomial  for  the  in¬ 
terval  How ever,  tble  polynomial  will  oontatn  terns  of 

the.  29th  degree  and  thus  14  operations  of  addition  sod  1$  opera¬ 
tions  of  aultlplloatlon  will  be  required  for  the  oaloulatlons. 
Proa  the  point  of  view  of  the  rate  of  oaloulatloa  of 
this  function  for  ^  £|xj  <  1,  it  is  more  advantageous  to  use  the 

following  formulas 

*  ! 

X 

arcais*  -  arctg y====— -  ,  (54) 


1 


from  whloh  the  oalculation  of  the  values  of  arcslnx  requires 
the  performance  of  8  addition,  7  multiplication,  2  division  and 
3  comparison  operations. 

Although  formula  (54)  holds  for  all  |x|<l,  for  (xi  <  ^ 
It  la  less  useful  than  formulas  (51)»  (52)  and  (53). 


7.  Calculation  of  the  function  f(x)  »  ^ 

Some  computers  do  not  have  the  division  operation; 
hence  an  Iteration  prooess  Is  used  to  obtain  Inverse  quantities. 

We  shall  examine  the  possibility  of  obtaining  Inverse 
quantities  using  best  polynomials  constructed  for  different 
segments  Into  whloh  the  Sntlre  Interval  of  definition  of  the 
Independent  variable  Is  divided. 

Ve  shall  assume  that  the  Inverse  quantity  Is  calcula¬ 
ted  for  a  normalised  argument,  l.e.  that 

We  shall  divide  the  intire  Interval  of  definition  of 
the  Independent  variable  Into  32  partial  segments,  and  at  the 
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midpoint  of  eaoh  of  those  segments  we  shall  calculate  the  val- 
uea  of  the  funotlon  and  Its  first  6  derivatives. 

Through  estimations  we  establish  that  funotlon  f(x)sl. 
oan  be  approximated  aoeurately  to  10  deolmai  plaoes  by  fifth- 
degree  polynomials  of  the  form  (13)  on  the  16  left  segments,  and 
by  fourth-degree  polynomials  of  the  fora  (12)  on  the  16  right 
Intervals.  The  coefficients  of  the  polynomials  are  oaloulated 
just  as  they  were  oaloulated  for  polynomials  (12)  and  (13). 
Calculation  of  f(x)  =  ^  with  these  polynomials  will  require  6 
addition  and  5  multiplication,  or  5  addition  and  4  multiplica¬ 
tion  operations.  Additionally,  5  conditional  branohes  will  be 
necessary  to  search  the  required  Interval,  in  the  computer's 
memory  unit  208  constants  must  be  stored. 

In  oonolualon  we  shall  determine  the  values  A^  and  Ag 
for  the  methods  examined  above  of  calculating  the  elementary 
functions.  This  will  permit  an  estimation  of  the  lncrese  In  the 
rate  of  calculation  of  these  functions.  In  oases  where  the  de¬ 
gree  of  the  approximating  polynomial  does  not  lnorease  from  one 
sub-interval  to  the  next  (In  calculating  ex  by  the  second  method, 
^x  ,  arotgx),  the  value  of  A ^  will  be  determined  relative  to  the 
old  formulas  for  calculating  elementary  functions. 

In  finding  Ag  It  should  be  noted  thatTQ  (the  time  for 
searching  In  the  table)  usually  does  not  exceed  two  X  of  addi¬ 
tion,  since  In  the  majority  of  cases  we  have  a  division  into 
four  sub-intervals,  and  thus  two  comparisons  must  be  performed 
to  find  the  necessary  Interval,  in  the  ease  where  the  number  of 
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sub- Intervals  la  greater  than  four,  the  nuaber  of  ooaparlaona 
will  be  ln2i,  where  t  la  the  nuaber  of  sub-lntervalsfTJ,  Typi- 
oally,  In  all  oaaaa  the  boundaries  of  the  sub-intervals  can  be 
defined  bj  two  or  three  deolaal  plaoea.  This  will  permit  the 
lnaertlon  Into  the  computer  oode  of  an  lnatrhotldnr for'  shortened 
comparison  whloh  oan  effect  a  comparison  over  half  of  a  word  In 
the  oonputer.  Additionally,  If  autoaatlo  aearoh  of  the  table  la 
employed,  as  in  the  IBM-650  oonputer  [ll,  or  the  arithmetic  unit 
la  specially  adapted  for  suoh  problems,  the  searoh  time  oan  be 
out  at  least  In  half.  Henee  we  shall  estimate  A2  from  two 
sides.  The  right  estimate  will  be  suitable  for  computers  not 
specially  adapted  for  the  problem  of  searohlng  In  the  table, 
while  the  left  estimate  will  be  suitable  for  computers  whloh 
are  thus  adapted.  Here  all  estimations  will  be  based  on  the 
assumption  that  the  ratio  of  the  time  required  to  perform  a 
multiplication  to  the  time  required  to  perform  an  addition 
operation  will  equal  two,  and  the  ratio  of  the  time  required 
to  perform  a  division  operation  to  the  time  required  to  perform 
an  addition  operation  will  equal  four.  These  ratios  hold  In  i 
best  present-day  computers.  If  these  ratios  are  larger,  the 
values  of  A^  and  A2  will  diminish  significantly. 

I 

The  application  of  subroutines  whloh  use  the  methods 
discussed  above  to  oaloulate  elementary  functions  requires  an 
additional  volume  of  memory  units.  This  lnorease  In  volume  Is 
oaused  by: 

1)  the  additional  storage  of  constants  for  polynomials 
whloh  approximate  the  value  of  a  funotlon  In  various  sub-lnter- 
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vals,  and  of  constants  which  define  the  boundaries  ef  these  aub- 
lntervalej 

2)  lnatruotlona  necessary  for  tba  oaleulatloa  of  tha 
polynomials  In  thaaa  sub-intervals; 

3)  lnatruotlona  which  affect  a  branoh  to  a  sub- routine 
for  calculating  tba  polynomial  in  tha  appropriate  sub-interval 
(in  the  simplest  oaae  these  will  be  instructions  for  comparing 
values  of  the  argument  with  the  table  of  sub-lntefval  boundaries). 

An  analysis  of  the  sub-routines  in  the  BESM-1  code  indi¬ 
cates  that  the  sub-routines  are  Increased  as  follows) 


C-alniUtird 

function 

BSEM-1  nib-routine*  { 

New  tub-routine 

No. 

No.  of 
con- 
ttanti 

■ 

Cell* 

occu¬ 

pied 

No.  of 
coo- 
itantt 

No.  of 
in*t« 

Cell* 

occu¬ 

pied 

1 

si  rue,  cos  x, 
tj'JC 

12 

38 

SO 

33 

40 

73 

2 

lnx 

8 

19 

27 

28 

22 

50 

3 

0*  lit  method 

14 

25 

39 

33 

28 

61 

g*  2nd  method 

- 

- 

- 

.96 

21 

117 

\/x  lit  method 

10 

27 

37 

102 

28 

130 

4 

\Jx  2nd  method 

- 

- 

- 

54 

24 

88 

\Jx  3rd  method 

- 

- 

- 

80 

25 

105 

5 

arctgx 

13 

23 

41 

18 

25 

43 

6 

arcsiiuc 

24 

57 

81 

27 

67-  ; 

94 

7 

10—2 

12 

25 

37 

12 

25 

37 

8 

2—10 

13 

26 

39 

13 

26 

39 

Total 

S! 

351  | 

ft 

540 

-54- 


CONCLUSIONS 


Ve  hart  desorlbed  a  net  bod  which,  at  the  expense  of  a 
oertaln  additional  loading  of  the  computer's  memory  unit,  per- 
nlta  an  inorease  in  the  rate  of  eeloulation  of  elementary  funo- 

r  * 

tlona. 


The  formulae  suggested  here  are  not  the  best  for  every 

.  *»  ^ 

eleotrorilo  ooaputer.  For  any  speoifie  ooaputer,  appropriate 
formulas  should  be  se looted  whloh  allow  for  the  properties  of 
the  ooaputer  relative  to  the  tine  required  for  the  performance 
of  arithmetic  operations,  whloh  Insure  oaloulatlon  with  the 
maximum  speed,  and  whose  programming  does  not  severely  load  the 
computer's  memory  unit. 

The  formulas  presented  here  can  be  quite  suitable  for 
specialized  computers,  where  the  occurrence  of  values  of  the 
function's  argument  -la  not  uniform  in  the  reduoed  Interval. 

In  this  case.  Instead  of  using  all  of  the  divisions,  we  may 

- 

use  only  certain  ones  and  thus  lnorease  the  rate  of  operation 
of  the.  oomputer  with  only  slight  additional  loading  of  the 
memory  unit. 

If  the  methods  discussed  in  3  and  6  are  used,  the 
number  of  multiplications  in  the  formulas  presented  here  can 
be  reduoed  still  further  because  of  the  increase  in  the  number 
addition  operations. 


All  formulas  derived  here  are  designed  for  calculation 
of  functions  aoeurata  to  10  decimal  places. 

If  th«  value a  of  funotlons  ara  to  be  calculated  with 
fewer  digits,  the  formulae  may  be  simplified  considerably. 

One  of  the  Indicated  variants  for  calculating  elemen¬ 
tary  functions  will  require  tbe  use  of  approximately  200  memory 
units  to  store  necessary  constants.  Tbe  maximum  number  of 
constants  required  for  the  aocompliabment  of  all  of  the  calcu¬ 
lation  sub-routines  for  tbe  elementary  functions  examined  here 
Is  around  300.  However,  around  200  of  these  are  used  to  formu- 
late  the  sub- routine a  of  the  funotlon  V?  with  4th-  and  5th- 
degree  polynomials,  and  the  function  ex  using  4th-degree  poly¬ 
nomials. 


The  use  of  the  suggested  formulas  for  the  fbnotlon 
arotgx  significantly  increases  tbe  rate  of  calculation  of  the 

I 

arctangent. 

The  methods  desorlbed  here  apparently  can  serve  as  one 
way  of  achieving  a  significant  reduction  of  the  operating  time 
of  the  computer  during  the  solution  of  various  problems. 

These  methods  can  be  transferred  In  their  entirety  to 
the  case  of  rational  approximation  of  elementary  functions.  In 

this  oase  their  application  permits  a  significant  eoonomy  of 

►  * 

computer  operating  time. 
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1.  Throughout  the  article  we  use  the  phrase  "beat 
polynomials  to  denote  poljnoalala  obtained  from  the  Taylor 
expansion  of  a  function  bj  multiple  rotation  with  Chebysher 
'*  polynomials.  These  polynomials  are  dose  to  polynomials 
which  are  best  in  the  generally  aooepted  sense.  The  error 
in  approximating  polynomials  used  in  the  article  conslats  of 
the  remainder  term  of  the  Taylor  series  and  the  errors  whioh 
arise  with  eaoh  rotation. 
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